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We derive the phase space density of bosons from a general boson interfer- 



two-particles and the single particle density distribution functions. If the bo- 
\^ ' son density is large, the two particles density distribution function can not 
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ever, if the boson density is so small that two particles density distribution 
function can be expressed as a product of two single particle density distribu- 
tions, then Bertsch's formula is recovered. For a Gaussian model, the effects 
of multi-particles Bose-Einstein correlations on the mean phase space density 
are studied. 
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I. INTRODUCTION 

The principal aim of the study of relativistic heavy-ion collisions is the search for evi- 
dences of the state of a quark gluon plasma (QGP) in the early stage of the reactions [1,2]. 
A quantity of great interest for the study of the QGP is maximum energy density that has 
been reached in the experiment. This maximum energy density is connected with the final 
state energy density and phase space volume. For ultra high energy collisions, more than 
eighty percent of final state particles are pions; therefore, it is very important to measure 
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the phase space density of pions which can be used to infer the energy density in the early 
stage of the coUisions [3]. Besides this, it is very important to estimate the phase space 
density also for the following reasons: (1) If the boson phase space density is very large, 
then pions will tend to stay in the same state and pion condensate may occur. (2): if the 
density of pions is very large, then the mean free-length among pions will be small; therefore, 
the evolution of pions in the final state should be described by the hydrodynamical equation 
[4]. 

Bertsch suggested a method which uses pion interferometry measurements [5] to calculate 
the mean pion phase space density several years ago [6]. Since then several calculations have 
been done at AGS energy [7] and SPS energy [6,8-12]. It has been found that the phase 
space density of pions is very low at AGS and SPS energies [6-8,10,12]; however this may 
be not the case at RHIC energies. 

Bose-Einstein (BE) correlation effects on the pion multiphcity distribution, on the single 
pion distribution, and on two-pion interferometry have been studied by many authors for a 
Gaussian source distribution [13-17]. It has been shown by Bialas and Zalewski that those 
results are valid for a wide class of models [18-20] . Fialkowski and Wit have implemented 
multi-particle Bose-Einstein correlations in Monte-Carlo generators and have studied Bose- 
Einstein correlation effects on the W mass shift, pion multiplicity distribution [21,22]. But 
the effects of multi-particle Bose-Einstein correlations on the mean phase space density have 
never been studied before. Theoretical studies [17,16] have shown that pion interferometry 
depends strongly on the pion multiplicity distribution which was overlooked in previous HBT 
analyses. Thus it is interesting for us to study the effects of the general pion interferometry 
formula on the mean phase space density. 

This paper is arranged in the following way. In the Section II, we re-derive Bertsch's 
formula and point out its implicit assumption. In Sec. Ill, we derive a phase space density 
formula from the general pion interferometry and we find that if the phase space density 
becomes small this new expression will be the same as the Bertsch's except a extra normal- 
ization factor. Unfortunately this simple relationship does not hold when the phase space 
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density becomes large. In Sec. IV, multi-pion Bose-Einstein correlation effects on the mean 
phase space density are studied. We find that multi-particle BE correlations will increase 
the mean phase space density. Finally we give our conclusions in Sec. V. 

II. BERTSCH'S FORMULA AND IMPLICIT ASSUMPTIONS 

The two-pion interferometry formula can be written as 



Here k — (pi + p2)/2 and q — pi — P2 are two-pions average momentum and relative 
momentum respectively. g^{x, k) is a Wigner function which can be interpreted as the 
probability of finding a pion at position x with momentum k. P2(Pi,P2) and Pi{p) are 
two-particle and single-particle inclusive distributions which are defined by 
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From Eqs. (1, 2), it is easily checked that if [6,10] 
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and 



P2(P1,P2) - Pl{pi)Pl{p2) = 7^ / d'xd'yf{x,k) 



{2ny 

f\y,k)exp{iq{x - y). (6) 
The reason that we put ^2:^^- in Eq. (4) is that in statistical physics for a infinite volume [24] 



Here V is the volume and f{p) is the Bose- Einstein distribution. 
Integrate Eq.(6) over q we have 



(27r)= 

The average phase space density {f^)k can be defined as 

, _ Jf{x,kfd'x 
Jf^{x,k)d^x 
_ J d^q[P2{k +I,k-I)-Pi{k + l)Pi{k - 1; 



Pi{k) 

Jd'q[Ci{q,k)-l]P,{k+l)P,{k-l) 



(9) 



Pi(k) 

Using the smooth approximation, pi ~ P2 ~ k, which has been shown to be valid in heavy-ion 
collisions for its large phase space [25], we have 

{f)k^Pi{k) J d'q[C^ik,q)-l]. (10) 

In Refs. [7,8,10], the authors have calculated the phase space density by assuming that 

Ci{pi,P2)^l + 

Xexp{~qlRl - ]^qlRl - ^qfPf - 2RoRiqoqi). (11) 

But the above parameterization of two-pion interferometry (Eqs. (11)) is not general, as 
in practice the two-pion correlation is fitted using function [16,15] 
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Cr{Pl,P2)^AC^ip^,P2). (12) 

Here ^4 is a normalization factor which exists in the two-pion interferometry formula [16,15]. 
If we use C|^(q', k) to take the place of C|(g, k) in Eq. (10), the phase space density will be 

{f)T^{A-l)P,{k)+Ax{f%. (13) 

Thus the phase space density will increase if A is larger than one or decrease if A is smaller 
than one. In the latter part of this paper we will show that this extra factor A though has 
been used in the data analyses will not appear in the phase space density formula on the 
condition that the phase space is large and the density is small. This guarantee that the 
application of Bertsch formula for heavy-ion collisions system is appropriate if the phase 
space density is small. 

III. PHASE SPACE DENSITY FROM THE GENERAL PION 
INTERFEROMETRY FORMULA 

It has been shown in Refs. [17,16] that the single particle spectrum, the two-particles 
spectrum and the two-pion interferometry formula read [23] 

^ max 

Pi{p) = E h^G,{p,p), (14) 

i=l 

^max 1 ^max ^ 

P2{Pl:P2) = J2 J2 hi+j[Gi{pi,Pi)Gi{p2,P2) 

i=l j=l 

+Gi{pi,P2)Gi{p2,Pi)]. (15) 

Where the definitions of hi and Gi{p,q) can be found in Ref. [17]. N^ax is the maximum 
multiplicity in the experiment. If N^ax — oo, we will obtain the formula of ^2(^1,^2) and 
Pi{p) given in Refs. [17,16]. 

The two-pion correlation function is [17] 

^2(^1,^2) 



C2{Pl:P2) 



Pl{pi)Pl{p2) 

Cr(Pl,P2)[l + i?2(Pl,P2)] (16) 



with 



ES=T/ii^.Gi(pi,Pi)G,(p2,P2) ' ^ ^ 

and 

R2(P1,P2) = 

E£— Ef—-^ h,^^G,{p,,p,)G^{p,,p,) 

Ef=r-'Ef=r-'h,+jG,ip,,p,)Gjip2,P2y ^ ^ 

In Ref. [17], we have shown that R2{k, q)\q=oo — 0. So if the two-pion correlation function is 
expressed as Eq. (12), then 

A=cr(Pi,P2), (19) 

which is a function of q and k. We can always define a Wigner function 5'(a;, k) which fulfils 
the following equation [14,15,17] 

^ max f 

^ hiGi{pi,p2) = / S{x,k)exp{iqx)d'^x. (20) 
i=i 

Thus, 

Pi{p) = E /^*G^«(P,P) = / S{x,p)d^x. (21) 

i=l 

It has been shown in Refs. [17,14,15] that for a special multiplicity distribution p„ = 

Lo(n) 



^n^a^' , , or for a small phase space density and Pn is a Poisson form, we have hi^j = hihj. 
Then Eqs. (15, 16) change to (for Nmax = oo) [14,15,17] 

P2(Pl,P2)-Pl(Pl)Pl(P2) + 

J S{x, k)S{y, k) exp{iq{x - y))dxdy, (22) 



and 



C2(Pl,P2) = 1 + 

/ S{x, k)S{y, k) exp{iq{x — y))dxdy 
JdxdyS{x,pi)S{y,P2) 



(23) 
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So we obtain Eq. (1). Thus all the derivations given in Ref. [6] are valid. This implies that 
the normalization factor A must be one and this can be shown from Eq. (17) under the 
condition that Nmax = oo and hi+j = hihj. From Eq.(13), we have 

{f)T - {f)k. (24) 

This verifies the rightness of the application of Bertsch's formula for the case of small boson 
densities. However, the relationship hi+j — hihj does not hold for all cases. For the following 
four kinds of multiplicity distributions: 

Vn = — j-exp(-(n)) 

77/. 

{Poisson distribution), 

(n)" 



Pn 



Pn 



(1 + (n))"+i 

{Bose — Einstein distribution), 

(n + fc-1)! ((n)/fc)" 
n\{k-l)\ (i + M)n+fc 

{negative binomial distribution), 

{Gamma distribution), (25) 

We can prove that [17] that if the phase space volume is large, h2/h\ ~ h^/{hih2) ~ 1 for 
the Poisson distribution, h2/h\ — 2 and hs/{hih2) — 3 for the Bose- Einstein distribution. 

If there are strong BE correlations among bosons, then it is impossible to express the 
two-particles distribution S{x,pi,y,p2) as S{x,pi)S{y,p2)- However, one can find a real 
function Si{x, k) which fulfils the following equation 

Gi{p\,P2) ^ j Si{x,k)exTp{iqx)d'^x. (26) 

Then we define S{x, k; y, k) and S{x,pi; y,P2) as 

S{x,k;y,k)^ £ £ h+jSi{x,k)Sj{y,k), (27) 



and 



S{x,pi;y,p2) 



^max 1 ^max i 

E E 

hi+jSi{x,pi)Sj{y,p2), 



(28) 



which satisfies 



^max 1 ^max 



E E hi+jGi{pi,p2)Gj{p2,Pi) 

/;;,.;:.)e.«.-.))... 



(29) 



and 



E E hi+jGi{pi,pi)Gj{p2,P2) 



i=l 



y" S{x,pi;y,p2)dxdy. 



(30) 



Because p*{pi,P2) = p{P2,Pi), then G*{p,p) = Gi{p,p) and G*{p,q) = Gi{q,p). Thus 
S{x,pi;y,p2) must be a real function which fulfils the requirement of the Wigner function. 
If /li+j = hihj and A^maa: — CO, from Eq. (28), Eq. (26) and Eq. (20), we have S{x,pi; y,p2) — 
S{x,pi)S{y,p2), thus Eqs. (21,22,23) are obtained. So we can identify S{x,pi,y,p2) as a 
two pion distribution function which represents the probability of two pions emitted from 
point X with momentum pi and from point y with momentum p2 respectively. Two-particle 
spectrum distributions can be written as 



-P2(pi,P2) = J S{x,pi;y,p2)dxdy + 

j S{x, k; y, k) exp [iq{x — y)]dxdy. 



(31) 



Then i?2(pi,P2) reads 



R2{Pl,P2) = 



JS{x, k; y, k) exp{iq{x - y))dxdy 
J S{x,pi;y,p2)dxdy 



(32) 



If [6] 



S{x, k; y, k) 5{xo - T)6{yo - T)f{x, k; y, k) 



(33) 



and 



S{x,k) ^ 5{xo-T)-^J{x,k), (34) 



then 



(35) 



. . ^ J f{x, k; X, k)d'^x 
Jf{x,k)d^x 

Using the smoothness approximation, Pi(pi) ~ P\{p2) ^ Pi{k), we have 

= A{f%. (36) 

In the above we have used the approximation A ~ C2^^{q,k) ~ const. It is interesting to 
point out that normally A should be a function of q and A;; however, in heavy-ion collisions, 
the practice is to normally fit it as a constant if the phase space density is small [17]. If the 
phase space density is high, it has been suggested in Ref. [17] to fit data using the function 

C2(P1,P2) = Cr{Pl,P2)[l + R2{P1,P2)]. (37) 

Here C^^'{q,k) = M[l + B{k) ■ exp{-q^Rl^^{k)) and R2 = X{k) exp{-q'^R'^{k)). Thus it is 
not the best choice to use a constant A in Eq.(12). 

IV. MULTI-PARTICLE BE CORRELATION EFFECTS ON THE MEAN PHASE 

SPACE DENSITY 

In the following we will study the effects of multi-pion BE correlations on the mean phase 
space density. We assume g\x,p) to be [13,26] 

2 2 
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then f^{x,p) reads 

Due to Eq. (5), we immediately come to the conclusion that no is the mean pion multiplicity 
observed in the experiment. It is easily checked that 

However, we have neglected the high-order BE correlation effects to get this {f^)k- H we keep 
only the leading terms in Eqs. (19,20) (correspondingly, this implies that we have assumed 
that the phase space volume is large), then 

Pi(p) = /iiGi(p,p), (41) 

and 

P2{Pl,P2) = h2[Gi{pi,pi)Gi{p2,P2) + 

Gi{pi,P2)G,ip2,Pi)]. (42) 

Thus the corresponding two-pion correlation function becomes 

^/ N ^2r, , Gi{pi,p2)Gi{p2,Pi) ^ , . 

C2{Pl:P2) = T2 1 + TTT VTTl T • (43) 

hi Giipi,pi)Gi{p2,P2) 
Bring Eq. (43) into Eq. (9) and ignore the extra normahzation factor, we get Eq. (40); 
furthermore, if we bring Eq. (43) into Eq. (35) and ignore the extra normalization factor, 
Eq. (40) will be regained. This results demonstrate that the Bertsch formula will give the 
correct result on the condition that the phase space volume is large. In the following we 
will study the multi-particles correlation effects on the mean phase space density due to the 
large phase space density. 
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FIG. 1. ri{k) as a function of k/A. Here the mean pion multiplicity is 20. The solid line, 
dashed line, dotted line and dot-dashed lines correspond to Poisson, Gamma, Negative binomial, 
and Bose-Einstein distribution. The phase space v = 2RA = 4 and 16 respectively. 



We define a function ri{k) as 



ri{k) 



(44) 



In the Fig.l, ri{k) is shown as a function of k/A. It is clear that for a large phase 
space, ri{k) ~ constant; on the other hand, if the phase space is small, {f)k and {f^)k 
have big differences at small momentum and small differences at large momentum. This 
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is easily understood since quantum effects are big for small momentum particles. It is 
interesting to notice that when the phase space is large, ri ~ /i2/^i(~ C^'^^iq, k)), which is 
two for the Bose-Einstein distribution, one for the Poisson distribution. From Eq. (33) and 
Eq. (29) (taking g = 0), we have 



/(x, k; y, k)d xd y = (27r) / S{x, k; y, k)d xd y 



(27r)6 E E K^,G,{k,k)G,{k,k). 
i=i j=i 



From Eq. (34) and Eq. (20) (taking g = 0), we get 



/(x, k)dx = {27Tf J S{x, k)d^x 

^max 

= {2Tif E hiG,{k,k). 
1=1 



From the definition of C2^'^{pi,p2)(Eq. (17)) and Eqs. (45, 46), we have 

J f{x,k;y,k)d^xd^y 



criq.k),=o 



J d^xd^yf{x,k)f{y,k)' 



(45) 



(46) 



(47) 
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FIG. 2. C2^^{q,k)\q=Q vs. k/A. Here the mean pion multiplicity is 20. The solid line (circles), 
dashed line (squares), dotted line (diamonds) and dot-dashed lines (triangles) corresponds to Pois- 
son, Gamma, Negative binomial, and Bose-Einstein distribution. The phase space v = 2RA = 4 
and 16 respectively. 
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In Fig. 2, Cl'^^{q,k)q=Q vs A;/ A is shown. We notice that when phase space is large, 
C2^^{q, k)q=o is a constant. In this case, / d^xd^yf{x, k;y,k) = 2 J d^d^yf{x, k)f{y, k) for the 
Bose-Einstein distribution and / d^xd^yf{x,k;y,k) = J d^d^yf{x,k)f{y,k) for the Poisson 
multiphcity distribution. But these relationships do not hold anymore for a small phase 
space volume. 

In Eqs. (9,36), we find that {f)k is A times larger than {f^)k when the phase space 
density is small and the function form of {f)k and {f^)k are different in the numerator. In 
the following we would like to show the relationship between f{x, k; x, k) and f{x, k). It is 
found that 



(2 



J fix,k;y,k)e'^^--y^dVxd'y 
{27rf J S{x, k] y, k)e"'^''-y^d^xd*yd^q 

E E 

1=1 j=i 

hi+jGi{pi,p2)Gj{p2,pi) (48) 



and 



f{x,k)d'x = i27rr d\ E E 

i=i j=i 

hihjGi{pi,p2)Gj{p2,Pi). (49) 

In the above derivation, we have used Eqs. (20, 26, 29, 33, 34). We define a function r2{k) 
as 

/7\ J f{x, k; x, k)d'^x , . 

''^^^ ^ Jd^xPix^k) • 

Similar to Ref. [17], we can prove that 



r2{k) ~ C2^^{q, ^) ~ r| = constant v — >• oo . (51) 
hi 



In the following, we will discuss the effects of multi-boson correlations on the distribution 
function of pions in the momentum space. If the multi-boson symmetrization effects are 
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small, the distribution function is f^{x,p). If the multi-boson correlations are strong, the 
distribution function is f{x, k). We define a function 

' Jf'{x,k)d^x Jf'{x,k)d^x' ^ ' 

which reflects the effects of multi-boson correlations on the distribution function. In Fig.3, 

r^ik) vs. /c/A is shown. It is interesting to notice that when phase space is large, 

r^{k) = — > cxD (53) 

for all distributions. On the other hand, when phase space is small, boson density becomes 
large at small momentum region. We define {f)l. as 

I Jf{x,k)dx 

This definition is similar to Eq. (9) but with f{x,k) taking the place of f^{x,k). From 
Eq. (9), Eq. (12), Eq. (35), Eq. (36), and Eq. (51), we get 



{f)k 

r2{k) 

C2{q,K) 



ik) J d^[^^^ -l]v^oo. (55) 



This is one of the main results of this paper. This result guarantees that the application 
of Bertsch's formula for a large system and dilute gas is appropriate. On the other hand, 
Bertsch formula is incomplete when the phase space density is high. This comes from the 
fact that we can not calculate / p{x^k)d^x from two-pion interferometry formula though 
we could calculate / /(x, k)d^x from the single particle spectrum. If the phase space density 
is high, C2^*(pi,p2) is not a constant anymore, we can not find the approximation formula 
as Eq. (55); however Eq. (35) can still be used to find the ratio of the number of particles 
pairs which are emitted from the same phase space cell to the average number of particles. 
The relation among (/)fc, and (/^)fc will be very complex when the phase space density 
is large. Then the physical meaning of the Bertsch formula is no longer clear. Multi-pion 
BE symmetrization will affect the current formalism and we believe that those effects will 
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be similar to the effects of multi-pion BE on the Wigner function g{x,p), which have been 
presented in Ref. [15] (or f^{x,p) presented here). 
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FIG. 3. r3(/c) as a function of k/A. Here the mean pion multiplicity is 20. The solid line 
(circles), dashed line (squares), dotted line (diamonds) and dot-dashed lines (triangles) corre- 
sponds to Poisson, Gamma, Negative binomial, and Bose-Einstein distribution. The phase space 
V = 2RA = 4 and 16 respectively. 



V. CONCLUSIONS 

In this paper, the mean phase space density distribution of bosons is derived from the 
general pion interferometry formula. We find that when the phase space is small and thus the 
boson density is high, the two particle source distribution can not be expressed as a product 
of two single particles source distributions. On the other hand, when the phase space is 
large and thus the boson density is small, Bertsch's formula is recovered. Thus Bertsch's 
formula can be used for the heavy-ion system if the freezeout pion phase space density is 
small. Multi-pion BE correlation effects on the mean phase space density distribution are 
studied, it is found that when the phase space density is large, bosons are concentrated in 
small momentum region and this effect is connected with the pion multiplicity distributions. 
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